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Abstract 



Starting from a simple marginally stable model considered for Lya- 
punov based boundary control of flexible mechanical systems, we add 
a term driving an instability and prove that for an appropriate control 
condition the system can become Lyapunov stable. A similar approx- 
imate extension is found for the general energy principle of linearized 
magnetohydro dynamics. The implementation of such external instan- 
taneous actions may however, impose challenging constraints for fu- 
sion plasmas. 
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1 Introduction 



Plasmas surrounded by resistive walls may be unstable within free boundary 
magnetohydrodynamic theory for fusion relevant /3 values (See [TJ [2]). To 
overcome this difficulty, a boundary control or feedback may be necessary 
(See [3]-[Tl]). Lyapunov methods have been applied successfully to the less 
demanding problem of boundary control of "flexible mechanical systems" as 
demonstrated in Refs.[T5]-[TT]. 

The simplest model of that kind [15] is governed by the equation of the 
displacement u(x,t) on the interval [0,L] of a flexible string: 

u tt ~u xx = (1) 

with u(0,t) = 0, u x (L,t) = —kut(L,t) with k > 0. To obtain a Lyapunov 
function for that model, we multiply Eq.(l) by u t and integrate over x from 
to L which leads, after integration by parts, to 

Sl^ + u ^ dx ~ [utUx] ° = ° - (2) 

Applying the boundary conditions above one obtains 

l§£ \\< + ul)dx = -ku t (L, tf. (3) 

The integral is positive and its time derivative is negative, which fulfills the 
condition for Lyapunov stability. However, the system is marginal for k = 
i.e. without boundary action. Since the relevant plasmas for higher /3 values 
are not only marginal but rather strongly unstable, the model (1) must be 
modified to permit an unstable situation in the absence of an external action 
which is the purpose of section 2. 



2 Unstable model 

The model now considered is a slight modification of Eq.(l) and is given by 

Utt - u xx - au = (4) 

with a > 0, u(0,t) = 0, u x (L,t) = —k\Ut{L,t) — k2u(L,t) where k\ and &2 
are positive constants. Multiplying Eq.(4) by u t we obtain 
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/ (ut + u 2 x — au 2 )dx — [u t u x }Q = 
Jo 



(5) 



l d r L 

2di 

after integration by parts over x. Inserting the new boundary conditions in 
(5) we have 



1 6 f L id 

2dtJo + U * ~ au ^ dx + ^dt^^ 1, ^ = ~ klUt ( L > *) ■ ( 6 ) 
We prove that 

k 



V= f [u 2 + u 2 x -au 2 + ^u 2 {L,t)]dx 
Jo L 



(7) 



is indefinite for fc 2 = if 7r/4 < L^fa < 37r/4 in the limit a = below, which 
by the second Lyapunov theorem means that the system is unstable without 
external action. Indeed, let us first minimize J L {u 2 — au 2 )dx for u(0,t) = 0. 
It leads to u xx + au = whose solution is u(x,t) = a(t)sin(xy/a). Then 
inserting the solution in (7) we obtain 

rL 

a sin 2 (x\^a) + a 2 acos 2 (x-Ja) — a 2 asin 2 (x\fa) 



Vrr 



■\~a 2 sin 2 {Ly/a) 

Li 



dx 



(8) 



or 



V„. 



d 2 sin 2 (x\fa) + c?acos{2x\fa) 
dx. 



+ ^-a 2 sin 2 (Ly/a) 

Li 



(9) 



Finally 



V m in 



L — 



sin(2Ly/a)\ a~\/(\ 



2^ J 
The system can be stabilized if 



+ 



sin {2Ly/a) + k 2 a 2 sin 2 {L^a). (10) 



ko > - 



a 
~2 



' L _ sm(2L v /o) ^ 



a 2 \/a 



2 v 'n J ' 2 

which in the limit of a = becomes 

k 2 > 



sin(2L\fa) 



(a 2 sin 2 (L\/a)j , 



a sin(2Ly/a) 



2 sin 2 (L^/a) 



(11) 
(12) 
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3 Magnet ohydro dynamic boundary control 



From the models above we have learned that k 2 is essential to influence 
positively the Lyapunov function. Let us now try to apply the procedure to 
the energy principle of linearized magnetohydrodynamics (See [U EJ [18]) 

^[(£•0 + SWp + SW P ,w + 5W w>inf ] = -a fjAfdS (13) 

where £ is the plasma displacement and the subscripts refer to the plasma, 
vacuum between plasma and wall and vacuum between wall and infinity 
respectively. At the plasma vacuum interface, we have nx A = (n-£)B where 
n is the normal to the interface, A is the perturbed vector potential at the 
vacuum side and B is the equilibrium magnetic field, a is the conductivity 
of the wall. The last volume integral in (13) can be transformed in a surface 
integral 

ISriSWwMf] = \t»I ( V x A ) 2dT = \ / n 5-AxVx AdS (14) 
2 at 2 at Jw.inf 2 Js 

where S is the external surface of the wall and is the normal to that 
surface. Eq.(13) can now be written 

^jrM-i) + SWp + 5Wp, w + [ n 5 • A x V x AdS] = -a [ (A) 2 dS (15) 

At this point it is usually assumed that x A = if the wall is taken as 
perfectly conducting. An alternative is to choose a boundary control law like 

n s xVxA— ^n 5 xVxA — k 2 A (16) 

on S with &2 constant in time but as a free function within S. Then (15) 
becomes 

n^ridO + SW P + 8Wp, w + 5W Wjinf + [ k 2 A 2 dS] = -a [ (A) 2 dS (17) 
2 at Js Jw 

Similarly to section 2 we can minimize, for £ = 0, 5W = 5Wp + 5Wp t w + 
5Ww,inf with respect to (£, A) with A vanishing for large distances to ob- 
tain (£ m m, A min ) and SW m i n . Since 5W m i n is indefinite the condition for 
stabilization is 

/ k 2 A 2 mm > SW mm . (18) 
Js 
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We can now determine k 2 ^ min , by optimizing the solutions of (18), as a func- 
tion of the position on the surface or, simply, assume a positive constant k 2 
satisfying (18). 

Let us notice that in the boundary law (16) a magnetic field is imposed 
on the surface S, which implies a surface current density. This is possible 
only if that surface is perfectly conducting. For resistive conductors the 
term J s k 2 A 2 dS will disappear in (17). Strictly speaking, a surface boundary 
control is not possible in the magnetohydrodynamic case unlike the case of 
section 2. In real situation the control on a boundary layer of (small) finite 
thickness seems necessary and such a treatment requires further investigation. 

4 Discussion 

An unstable model, in contrast to the marginal models of Ref.[15], can be 
stabilized by controlling the boundary if the external action (here k 2 ) is large 
enough. It turns out that the method used for this model can be extended 
to linearized magnetohydrodynamics for which codes exist to minimize the 
energy principle of Ref.[18] in the approximation that the external surface 
of the wall is covered by an infinitesimally thick layer of perfectly conduct- 
ing actuators. In addition to potentially extending the study by removing 
this assumption many problems remain open of how to realize instantaneous 
boundary conditions by coils in view of their induction, resistance etc. Start- 
ing from the sensors, the amplifyers should have high saturation levels and 
broad bands. Moreover, since the actuators must consist of coils, it means 
that they may have to be placed within the resistive vacuum vessel to avoid 
time delays. This may be a challenging problem for fusion reactors. 
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